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Abstract 

The production of leading neutrons, where the neutron carries a large fraction xl of the in- 
coming proton's longitudinal momentum, is studied in deep-inelastic positron-proton scat- 
tering at HERA. The data were taken with the HI detector in the years 2006 and 2007 
and correspond to an integrated luminosity of 122 pb _1 . The semi-inclusive cross section 
is measured in the phase space defined by the photon virtuality 6 < Q 2 < 100 GeV 2 , 
Bjorken scaling variable 1.5 • 10~ 4 < x < 3 • 10 -2 , longitudinal momentum fraction 
0.32 < xl < 0.95 and neutron transverse momentum pr < 0.2 GeV. The leading neutron 
structure function, Frf ^^(Q 2 , xl), and the fraction of deep-inelastic scattering events 
containing a leading neutron are studied as a function of Q 2 , x and xl- Assuming that 
the pion exchange mechanism dominates leading neutron production, the data provide con- 
straints on the shape of the pion structure function. 
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1 Introduction 



The production of leading baryons in deep-inelastic scattering (DIS) provides a testing ground 
for the theory of strong interactions in the soft regime. Events containing a neutron, which 
carries a large fraction xl of the longitudinal momentum of the incoming proton, have been 
observed in electron-proton collisions at HERA flU-H. The process leading to such events, 
ep — > e'nX, is illustrated in Fig. [T^. The pion exchange mechanism, illustrated in Fig. [Tt>, 
is expected to dominate leading neutron production at large x L and low transverse momentum 
of the neutron [[8UT41 . In this picture of leading neutron production, the proton fluctuates into 
a state consisting of a positively charged pion and a neutron p — > nn + . The virtual photon 
subsequently interacts with a parton from the pion. Consequently, the cross section factorises 
into two parts (proton vertex factorisation): one factor describes the proton fluctuation into 
a nix + state, the other describes the photon-pion scattering (Hill. The production of leading 
neutrons in DIS at HERA may therefore provide constraints on the structure of the pion at low 
to medium Bjorken-x, while the knowledge of the pion structure from fixed target experiments 
[fT5l is limited to higher x values. Previous HI and ZEUS studies of semi-inclusive leading 
neutron DIS cross sections 01121 demonstrate that these measurements are indeed sensitive to 
the structure of the pion and can distinguish between different parameterisations of the pion 
structure function. 

Energetic neutrons can also be produced in the fragmentation of the proton remnant. The 
comparison of leading neutron production with inclusive DIS provides tests of fragmentation 
mechanisms. The hypothesis of limiting fragmentation [fT6l[T7l states that, in the high-energy 
limit, the cross section for the inclusive production of particles in the target fragmentation re- 
gion becomes independent of the incident projectile energy. This hypothesis implies that, in 
DIS, leading neutron production is insensitive to Bjorken-a; and the virtuality of the exchanged 
photon Q 2 . 

In this paper a measurement of the semi-inclusive cross section for leading neutron produc- 
tion in DIS is presented. This analysis is based on a data sample corresponding to an integrated 
luminosity which is 36 times larger than that of the previous HI publication [Q]|. A new neutron 
calorimeter with improved performance is used. The larger data set together with better exper- 
imental capabilities allow the extension of the kinematic range of the measurement to higher 
values of Q 2 and x and the reduction of the total uncertainty of the results. 



2 Event Kinematics and Reconstruction 

The kinematics of semi-inclusive leading neutron production are shown in Fig. [Tk, where the 
four-vectors of the incoming and outgoing particles and of the exchanged virtual photon 7* are 
indicated. The kinematic variables Q 2 , x and y are used to describe the inclusive DIS scattering 
process. They are defined as 



Q 2 = -q 2 , 



Q 2 



p ■ q 

p ■ k 



(1) 



x = 



2p ■ q 



y = 



4 




where p, k and q are the four-momenta of the incident proton, the incident positron and the 
exchanged virtual photon, respectively. These variables are reconstructed using the technique 
introduced in [fT8l . which optimises the resolution throughout the measured y range by exploit- 
ing information from both the scattered positron and the hadronic final state. 

The kinematic variables used to describe the final state neutron are the longitudinal momen- 
tum fraction x L and the squared four-momentum transfer t between the incident proton and the 
final state neutron 

x L = l- q ' iP ~ Pn) ~ E n /E p , t={p-p n 
q-p 

where m v is the proton mass, p n is the four-momentum of the final state neutron, m n is the 
neutron mass and E n and pt are the neutron energy and transverse momentum, respectively. 

The four-fold differential cross section for leading neutron production can be parameterised 
by a semi-inclusive structure function, F^^, defined by 




d 4 a(ep — > enX) Ana 2 f ^ | y 
dx dQ 2 dx L dt ~ ~x~Q 4 



(l-y+ y -^)F 2 LNii \Q 2 ,x,x L ,t). (3) 



The contribution from longitudinally polarised photons can be neglected in the phase space 
studied in this analysis. Integrating Eq. [3] over t yields the semi-inclusive structure function 
F^ N ^ measured in this analysis 

d 3 a(ep ->■ enX) f tmin d A a(ep -> enX) 

~ Lit- 



dxdQ 2 dxL Ju dxdQ 2 dxLdt 



to 

4 ™ 2 ^- 3/ + fK~< 3 '(Q=, (4) 



xQ A V 2 



where the integration limits are 




^max\2 



i'unn ( 1 - •''/. ) ( — - mj t j and t = ^- L - + t min . (5) 



Here p™ ax is the upper limit of the neutron transverse momentum used for the F 2 mea- 
surement. Smaller values of p™ ax are expected to enhance the relative contribution of pion 
exchange HlOlllli In this analysis p™ ax is set to 0.2 GeV, as in the previous HI publication [T). 

3 Experimental Procedure and Data Analysis 

The data used in this analysis were collected with the HI detector at HERA in the years 2006 
and 2007 and correspond to an integrated luminosity of 122 pb^ 1 . During this period HERA 
collided positrons and protons with energies of E e = 27.6 GeV and E p = 920 GeV, respec- 
tively. 

3.1 HI detector 

A detailed description of the HI detector can be found elsewhere lfT9] - |23ll . Here, a brief account 
is given of the components most relevant to the present analysis. The origin of the right-handed 
HI coordinate system is the nominal ep interaction point. The direction of the proton beam 
defines the positive z-axis (forward direction); the polar angle 9 is measured with respect to this 
axis. Transverse momenta are measured in the x — y plane. 

The ep interaction region is surrounded by a two-layered silicon strip detector and two large 
concentric drift chambers. Using these detectors charged particle momenta are measured in the 
angular range 25° < 9 < 155° with a resolution of cr(p T )/p T = 0.005 p T /GeV © 0.015 Il24l . 
The tracking system is surrounded by a finely segmented Liquid Argon (LAr) calorimeter, 
which covers a range in polar angle of 4° < 9 < 154° with full azimuthal acceptance. The 
LAr calorimeter consists of an electromagnetic section with lead absorber and a hadronic sec- 
tion with steel absorber. The total depth of the LAr calorimeter ranges from 4.5 to 8 hadronic 
interaction lengths. Its energy resolution, determined in test beam measurements, is a(E)/E ~ 
12%/y£[GeV] © 1% for electrons and a(E)/E « 50%/y/E [GeV] © 2% for charged pi- 
ons ll2~5l . The backward region (153° < 9 < 177.8°) is covered by a lead/scintillating-fibre 
calorimeter, the SpaCal. Its main purpose is the detection of scattered positrons. The energy 
resolution for positrons is a(E)/E m 7.1%/ ^ E[GeV] © 1%. The LAr and SpaCal calorime- 
ters are surrounded by a superconducting solenoid which provides a uniform magnetic field of 
1.16 T along the beam direction. 

The luminosity is measured via the Bethe-Heitler process ep — >• e'pj. The final state photon 
is detected in a dedicated calorimeter situated near the beam pipe at z = —103 m. 

3.2 Detection of leading neutrons 

Leading neutrons are detected in the forward neutron calorimeter (FNC). The FNC is situated 
at a polar angle of 0° beyond the magnets used to deflect the proton beam, at z = 106 m. A 
schematic view of the FNC is shown in Fig. |2k. It consists of the Main Calorimeter and the 
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Preshower Calorimeter. In addition, two layers of veto counters situated at a distance of 2 m in 
front of the Preshower Calorimeter are used to veto charged particles. 

The Preshower Calorimeter is a 40 cm long lead-scintillator sandwich calorimeter, corre- 
sponding to about 60 radiation lengths or 1.6 hadronic interaction lengths. It is composed of 
24 planes: the first 12 planes each consist of a lead plate of 7.5 mm thickness and a scintil- 
lator plate of 2.6 mm thickness, the second 12 planes each consist of a lead plate of 14 mm 
thickness and a scintillator plate of 5.2 mm thickness. The transverse size of the scintillating 
plates is 26 x 26 cm 2 . Each scintillating plate has 45 parallel grooves holding 1.2 mm diameter 
wavelength shifter (WLS) fibres. In order to obtain good spatial resolution, the orientation of 
fibres alternates from horizontal to vertical in consecutive planes. At each plane the fibres are 
bundled into nine strips of five fibres. Longitudinally, the strips are combined into 9 vertical 
and 9 horizontal towers which are finally connected to 18 photomultipliers. 




Figure 2: (a) Schematic view of the FNC. (b) Layout of tiles on an active board of the Main 
Calorimeter; the position of readout fibres is indicated by dotted lines; the hatched area shows 
the geometrical acceptance window defined by the beam line elements. The position corre- 
sponding to 9 = 0° is also indicated. All dimensions are given in mm. 

The Main Calorimeter of the FNC is a sandwich-type calorimeter consisting of four identical 
sections. Each section is 51.5 cm long with transverse dimensions of 60 x 60 cm 2 and consists 
of 25 lead absorber plates of 14 mm thickness and 25 active boards with 3 mm thick scintillators. 
Each active board is made of 6 scintillating tiles with a transverse size of 20 x 20 cm 2 and 2 
tiles of transverse size 20 x 26 cm 2 , as shown in Fig. [2b- Each scintillating tile has a circular 
groove holding a 1 mm diameter WLS fibre which is attached via light connectors to two 1 mm 
diameter transparent fibres. Longitudinally, for each section the fibres from 25 tiles with the 
same transverse position are bundled into a tower and connected to one photomultiplier. There 
are then 8 towers in each section, making a total of 32 towers in the Main Calorimeter. The 
proton beam pipe is located in a rectangular space along the top of the calorimeter, as can 
be seen in Fig. [2l The total length of the Main Calorimeter is 206 cm, corresponding to 8.9 
hadronic interaction lengths. 
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All modules of the FNC were initially calibrated at CERN using electron beams with en- 
ergies between 120 and 230 GeV and hadron beams with energies between 120 and 350 GeV. 
The FNC was positioned on a movable platform which allowed the response of each module 
and tower to be measured separately. After this initial calibration, the FNC had an approxi- 
mately uniform response independent of impact position. The linearity of the energy response 
of the FNC was measured at HERA from beam-gas interactions during dedicated runs, when 
the proton beam was accelerated to five intermediate energies between 150 GeV and 920 GeV. 
The energy response of the FNC is linear to a precision of 3% and the hadronic energy resolu- 
tion is a(E)/E £3 Q3%/ ^/E[GeV] © 3%. The energy resolution for electromagnetic showers, 
which are always fully contained in the Preshower Calorimeter, is a(E) /E « 20% / y/ E [GeV]. 
For hadronic showers starting in the Main Calorimeter, the spatial resolution is a(x,y) ~ 
10 cm/ y/ E [GeV] © 0.6 cm. A better spatial resolution of about 2 mm is achieved for the 
electromagnetic showers and for those hadronic showers starting in the Preshower Calorimeter. 

The acceptance of the FNC is defined by the aperture of the HERA beam line magnets and 
is limited to neutron scattering angles of 6 n £0.8 mrad with an approximately 30% azimuthal 
coverage. The geometrical acceptance window of the FNC is indicated in Fig. [2j?. 

After the calorimeter was installed in the HI beam line, the stability of calibration constants 
was monitored using interactions between the proton beam and residual gas molecules in the 
beam pipe. The neutron energy spectrum was compared with the results of Monte Carlo sim- 
ulation based on pion exchange. From this monitoring, the time dependent variations of the 
calibration constants of order of a few percent were determined. Short term variations of pho- 
tomultiplier gain were monitored using LED signals. The LEDs were operated during empty 
bunches at a frequency of 0.8 Hz. The averaged LED signal responses were used to provide 
offline energy corrections applied during the reconstruction. 

The longitudinal segmentation of the FNC allows efficient discrimination of neutrons from 
charged particles or electromagnetic energy deposits. Charged particles are rejected using sig- 
nals from the veto counters. Events with an energy deposit only in the Preshower Calorimeter 
are identified as electromagnetic showers (photons, n° mesons). All other patterns of energy 
deposition are identified as hadronic showers, which are subsequently classified as either those 
starting in the Preshower Calorimeter or those starting in the Main Calorimeter, the latter having 
worse spatial resolution. In this analysis both hadronic shower types are used. 

3.3 Event selection 

The data sample was collected using a trigger which requires the scattered positron to be mea- 
sured in the SpaCal and at least 100 GeV energy to be deposited in the FNC. The trigger signal 
from the FNC is formed using an analogue sum of signals from the Preshower Calorimeter and 
the central towers of the Main Calorimeter. The trigger efficiency is above 97% in most of the 
analysis phase space, decreasing to 94% for x L < 0.4. 

The selection of DIS events is based on the identification of the scattered positron as the most 
energetic compact calorimetric deposit in the SpaCal with an energy E' e > 11 GeV and a polar 
angle 156° < 9' e < 175°. The energy weighted cluster radius is required to be less than 4 cm, 
as expected for an electromagnetic shower ll26ll . The ^-coordinate of the primary event vertex 
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is required to be within ±35 cm of the nominal position of the interaction point. The remaining 
clusters in the calorimeters and the charged tracks are combined to reconstruct the hadronic final 
state. To suppress events with initial state hard photon radiation, as well as events originating 
from non-ep interactions, the quantity E—p z , summed over all reconstructed particles including 
the positron, is required to lie between 35 GeV and 70 GeV. This quantity is expected to be 
twice the electron beam energy for DIS events without QED radiation. Furthermore, events are 
selected within the kinematic range 6 < Q 2 < 100 GeV 2 , 0.02 < y < 0.6 and 1.5 • 10~ 4 < x < 
3 • 10- 2 . 

Events containing a leading neutron are selected by requiring a hadronic cluster in the FNC 
with an energy above 275 GeV and a polar angle below 0.75 mrad. The cut on polar angle, 
defined by the geometrical acceptance of the FNC, restricts the neutron transverse momenta 
to the range pr < %l ■ 0.69 GeV. 

The final data sample contains 315960 events which satisfy these selection criteria. For the 
measurement of F^^, an additional requirement on the transverse momentum of the neutron 
Pt < 0.2 GeV is applied to enhance the relative contribution of pion exchange and to avoid 
having an x L dependent p T cut. The number of events selected with this additional requirement 
is 209150. 



3.4 Monte Carlo simulation and corrections to the data 

Monte Carlo simulations are used to correct the data for the effects of detector acceptance, 
inefficiencies and migrations between measurement intervals due to finite resolution and QED 
radiation. All generated events are passed through a GEANT3 ll2~7ll based simulation of the HI 
apparatus and are processed using the same reconstruction and analysis framework as is used 
for the data. 

The DJANGO ll28l program generates inclusive DIS events. It is based on leading order 
electroweak cross sections and takes QCD effects into account up to order a s . The hadronic final 
state is simulated using ARIADNE 11291 , based on the Colour Dipole Model, with subsequent 
hadronisation effects modelled using the Lund string fragmentation model as implemented in 
JETSET [30J. DJANGO is also used in this analysis to simulate events where leading neutrons 
originate from proton remnant fragmentation. RAPGAP OTI is a general purpose event gener- 
ator for inclusive and diffractive ep interactions. Higher order QCD effects are simulated using 
parton showers and the final state hadrons are obtained via Lund string fragmentation. Higher 
order electroweak processes in the DJANGO and RAPGAP generators are simulated using an 
interface to HERACLES 021 

In the version denoted below as RAPGAP-7T, the program simulates exclusively the scatter- 
ing of virtual or real photons off an exchanged pion. Here, the cross section for photon-proton 
scattering to the final state nX takes the form 

da(ep — > e'nX) = f n +/ p (xL,t) ■ dcr(e7r + — >■ e'X), (6) 

where /„-+ / p (xl, t) represents the pion flux associated with the splitting of a proton into a n + n 
system and da(e7r + — > e'X) is the cross section of the positron-pion interaction. There are 
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several parameterisations of the pion flux [|9]-[T3l. In this analysis, the pion flux factor is taken 
from the light-cone representation iTTOll as 

where is the pion mass, gl^/Ait = 13.6 is the pirn coupling constant deduced from phe- 
nomenological analyses of low-energy data [33] and R wn = 0.93 GeV -1 is the radius of the 
pion-proton Fock state IfTOl . 

The DJANGO and RAPGAP-7T Monte Carlo simulations are calculated using GRV leading 
order parton distributions for the proton 041 and the pion [35], respectively. 

Figure [3] shows the observed energy and p T distribution of the neutron for selected data 
sample together with Monte Carlo simulations. The differences between the RAPGAP-7T and 
DJANGO generators are particularly visible in the neutron energy distributions (Fig. [3k). The 
RAPGAP-7T simulation peaks at E n ~ 650 GeV and describes the shape of the distribution at 
high energies. At lower energies (E n ^600 GeV) the RAPGAP-7T simulation does not describe 
the data. In this region additional physics processes which are expected to contribute signifi- 
cantly are not simulated. In contrast, the DJANGO simulation predicts a large contribution to 
the cross section in this region. The best description of the data is achieved if the predictions 
of the RAPGAP-7T and DJANGO Monte Carlo programs are added, using weighting factors 
of 0.65 and 1.2 for RAPGAP-7T and DJANGO, respectively. This Monte Carlo combination is 
labelled as "0.65xRAPGAP-7r + 1.2xDJANGO" in the figures and is used to correct the data. 

Cross sections at hadron level are determined from selected events by applying bin depen- 
dent correction factors. These factors are determined from the combination of DJANGO and 
RAPGAP-7T Monte Carlo simulations as the ratios of the cross sections obtained from particles 
at hadron level without QED radiation to the cross section calculated using reconstructed parti- 
cles and including QED radiation effects. The typical value of these factors is about 1.2 at the 
highest Xl increasing to 4 at the lowest x L due to the non-uniform azimuthal acceptance of the 
FNC. 

The binning in Q 2 , x and xl used to measure is given in Tabled! The bin purities, 

defined as the fraction of events reconstructed in a particular bin that originated from that bin on 
hadron level, and the bin stabilities, defined as the fraction of events originating from a particular 
bin on hadron level that are reconstructed in that bin, are higher than 50% for (Q 2 , x)-bins and 
higher than 60% for x L bins. 

The measured distributions may contain background arising from different sources. The 
background from photoproduction processes, where the positron is scattered into the backward 
beam pipe and a particle from the hadronic final state fakes the positron signature in the SpaCal, 
is estimated using the PHOJET Monte Carlo generator [|36ll . This background is negligible ex- 
cept at the highest y values where it can reach 6%. Background also arises from the random 
coincidence of DIS events, causing activity in the central detector, with proton beam-gas inter- 
actions, which give a neutron signal in the FNC. This contribution is estimated by combining 
DIS events with neutrons originating from beam-gas interactions in the bunch-crossings adja- 
cent to the bunch-crossing of the DIS event. It is found to be smaller than 1%. The estimated 
background contributions are not subtracted from the measurements. 
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The contribution from proton dissociation, where the leading neutron originates from the 
decay of a higher mass state, is estimated using an implementation in RAPGAP of the dissocia- 
tion model originally developed for the DIFFVM [37J Monte Carlo generator. This contribution 
can be up to 30% at low xl values, but is less than 5% for xl > 0.7. It is included in the cross 
section definition. 



3.5 Systematic uncertainties 

The effects of various systematic uncertainties on the cross section measurements are deter- 
mined using Monte Carlo simulations, by propagating the corresponding estimated measure- 
ment uncertainty through the full analysis chain. 

The acceptance of the FNC is defined by the interaction point and the geometry of the 
HERA magnets and is determined using Monte Carlo simulations. The angular distribution of 
the neutrons studied in this analysis is sharply peaked in the forward direction. Therefore the 
acceptance depends critically on small inclinations of the incoming proton beam with respect to 
its nominal direction. The uncertainty on the position of the neutron impact point is estimated 
to be 5 mm, which results in a 4.4% uncertainty in the FNC acceptance. The uncertainty in the 
neutron detection efficiency and the uncertainty of 2% on the absolute energy scale of the FNC 
lead to a systematic error on the cross section of 2% and 6%, respectively. An additional 0.5% 
uncertainty is attributed to the trigger efficiency. These effects are strongly correlated between 
measurement intervals and mainly contribute to the overall normalisation uncertainty. 

The uncertainties on the measurements of the scattered positron energy (1%) and angle 
(1 mrad) in the SpaCal lead to a combined systematic uncertainty of typically 1.9% on the 
cross section. The uncertainty of the energy measurement of the hadronic final state in the 
central detectors affect the reconstruction of the kinematic variables y, Q 2 and x. This hadronic 
energy scale uncertainty is estimated to be 4% for this measurement, leading to an uncertainty 
on the cross section of 0.4% on average. 

The systematic error on the efficiency to reconstruct the event vertex is determined by com- 
paring the reconstruction efficiencies for the data and the Monte Carlo simulation. The discrep- 
ancy is less than 1%. 

The systematic uncertainty arising from the radiative corrections and the model dependence 
of the data correction are estimated by varying the DJANGO and RAPGAP- n scaling factors 
described in Sect. 13 .41 within values permitted by the data. The resulting uncertainty on the cross 
section is below 5% in most of the bins and typically 2%. 

The luminosity measurement uncertainty for the selected run period leads to an overall 
normalisation uncertainty of 5%. 

The total systematic error in each bin is calculated as the quadratic sum of all contributions. 
Systematic errors are typically 10% for F 2 iiV ^ and 14% for da/dxi measurements. 
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4 Results 



The single differential leading neutron DIS cross sections as a function of xl are presented in 
Fig. |4] and Table |2] In Fig. 0b the differential cross section in xl is shown for px < 0.2 GeV. 
The difference between the shapes of distributions in Fig. and Fig. @K in particular a steep 
fall at low xl in Fig. Hk, is due to the geometrical acceptance of the FNC which restricts the 
accessible p T range to p T < x L ■ 0.69 GeV. 

The measured cross sections are compared with the Monte Carlo simulations. For large 
values of xl £0.7, the RAPGAP-7T simulation describes the shape of the Xl distributions well, 
in agreement with the assumption that at high xl the dominant mechanism for leading neutron 
production is pion exchange. The full xl range is well described by the sum of the RAPGAP-7T 
and DJANGO Monte Carlo generators, using the scaling factors discussed in Sect. 13.41 This 
indicates that the n exchange mechanism dominates at high xl, while proton remnant fragmen- 
tation gives a significant contribution at low xl- 

The measured cross sections are also compared with the predictions of the Soft Colour In- 
teraction model (SCI) ll38l . implemented in the LEPTO Monte Carlo generator program ll3~9~l . 
In the SCI model, the production of leading baryons and diffraction-like configurations is en- 
hanced via non-perturbative colour rearrangements between the outgoing partons. A refined 
version of the model uses a generalised area law (GAL) 11401 for the colour rearrangement prob- 
ability. Compared to DJANGO MC predictions, SCI-GAL improves the description at higher 
xl, as can be seen in Fig. SI However, for lower pr values the predicted cross section for 
xl > 0.5 is still too low. 

Figure |5] and Tables [3l|6] present the measurement of the semi-inclusive structure function 
F 2 LN{3) (Q 2 ,x,x L ) in the range 6 < Q 2 < 100 GeV 2 , 1.5 • 10~ 4 < x < 3 ■ 10~ 2 , 0.32 < 
xl < 0.95 and pr < 0.2 GeV. In all (Q 2 ,x) bins, the shape of the F 2 LN ^ distribution as a 
function of xl is similar to the shape of the single differential cross section in xl (Fig.0)). The 
distributions are reasonably well described by the combination of RAPGAP-7T and DJANGO 
simulations. 

The measurement of F 2 LN ^ allows the validity of the hypothesis of limiting fragmentation 
to be tested, according to which the production of leading neutrons in the proton fragmentation 
region is independent of Q 2 and x. To investigate this prediction, the ratio of the semi-inclusive 
structure function F 2 LN to the proton structure function F 2 is studied as a function of Q 2 in 
bins of x and xl (Fig. [6l). The values of F 2 are obtained from the H1PDF2009 parameterisa- 
tion ll26ll . The horizontal lines indicate the average value of the ratio for a given xl bin. These 
average values decrease from 7% to 2% with increasing xl, reflecting the general behaviour 
observed in Figs. H] and [5l The ratios are almost independent of x and Q 2 in each xl bin, im- 
plying that F 2 LN and F 2 have a similar (Q 2 , x) behaviour, as expected from the hypothesis of 
limiting fragmentation. 

Assuming that leading neutrons are produced via the exchange of a colour singlet particle, 
e.g. the 7r + , the structure function F 2 LN ^ factorises into a flux factor which is a function of xl 
and a structure function F 2 LN ^ which depends on Q 2 and = x/(l — x£). The quantity /3 
can be interpreted as the fraction of the exchanged particle's momentum carried by the parton 
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interacting with the virtual photon. The value of F 2 LN<y?, \Q 2 , f3, xl) is obtained by replacing 
the variable x with (3 in Eq. HI The distribution of F 2 LN ^ 3 \ shown in Fig. [/J has a similar 
dependence on j3 in all (Q 2 ,xl) bins. This behaviour can be approximated by a power law 
function F 2 LN ^ oc f3~ x . In each (Q 2 ,xl) bin a fit of the parameter A is performed. Within 
uncertainties the value of the fitted parameter A is independent of xl which is consistent with 
proton vertex factorisation. 

However, as a function of Q 2 , the values of A increase from 0.23 at lowest Q 2 to 0.3 at 
highest Q 2 . This slow Q 2 dependence is similar to the rise towards low x of the proton structure 
function F 2 measured in inclusive DIS [1411 . It is further investigated by fitting the measured 
F% (Q 2 ) P, xl) assuming the following functional form 

F 2 LN{3) (Q 2 , 0, x L ) = c(x L ) • /T A(Q2) with A = a ■ HQ 2 /A 2 ) (8) 

where a, A and the normalisations c(xl) are the free parameters of the fit. This nine parameter 
fit over the whole xl range has a good x 2 , supporting the validity of the employed ansatz (Eq.[8]). 
Within the experimental uncertainties the obtained values of the parameters a = 0.052 ± 0.003 
and A = 0.416 ± 0.052 GeV are in reasonable agreement with those obtained in the analysis of 
the proton structure function F 2 (|4T|. This demonstrates the similarity between the Q 2 evolution 
of F 2 LN ^ and the Q 2 evolution of the proton structure function F 2 . 

Since pion exchange dominates leading neutron production at high xl and low pr, the mea- 
surement of F 2 LN ^ in the range 0.68 < xl < 0.77 can be used to estimate the pion structure 
function at low Bjorken-x, following the procedure introduced in [1J. Assuming proton vertex 
factorisation, which is supported by the present data as explained above, the quantity F 2 LN ^ /TV 
can be interpreted as being equal to the structure function of the pion, where 



T n (x L ) = / fw/ P (xL,t)dt (9) 



to 



is the integral of the pion flux over the measured t-range, where t and t min are given by 
Eqs. O The pion flux from Eq. [TJused for the RAPGAP-7T simulation yields T n = 0.13 for 
p™ ax = 0.2 GeV at xl = 0.73, which is the central value of the chosen xl range. Using other 
parameterisations of the pion flux, e.g. from ll9l [TT14T3l . leads to values of the pion flux integral 
which may differ by up to 30%. In this evaluation of the pion structure function contributions 
from background processes like the exchange of p and a 2 -mesons, proton diffractive dissocia- 
tion and A production are not taken into account. Within the narrow xl range considered here 
they are only expected to affect the absolute normalisation of the results. The contribution of 
neutrons from fragmentation is of order 25 — 35%, as can be estimated using the DJANGO 
prediction. The relative size of this contribution is largely independent of Q 2 and (3 and thus 
has little impact on the shape of the distribution. 

Figure [8] shows F 2 N /V v as a function of Q 2 in bins of f3 while Fig. [9] shows F 2 LN ^ /TV 
as a function of (3 in bins of Q 2 . In Fig. [9] the contribution of neutrons from fragmentation, 
as predicted by DJANGO and scaled by a weighting factor 1.2, as described in Sect. 13 .41 is 
indicated. The data are compared to predictions of the parameterisations of the pion struc- 
ture function GRSc-7r ll42l and ABFKW-7T [|43l . The measurements are also compared to the 
H1PDF2009 parameterisation of the proton structure function F 2 (Q 2 ,x) ll26l which is scaled 
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by a factor of 2/3 in order to naively account for the different number of valence quarks in the 
pion and proton, respectively. The values of F 2 (Q 2 , x) are calculated at Bjorken-x equal to f3. 
The Q 2 distribution exhibits a rise with increasing Q 2 (i.e. scaling violation) for all (3 values 
in the measured range, which is similar in size and shape to that seen in the parameterisations 
of the inclusive structure functions of both the pion and proton (Fig. [8]). The distributions 
show a steep rise with decreasing for all Q 2 values (Fig. [9]). This behaviour is in reasonable 
agreement with the pion and proton structure function parameterisations. In absolute values the 
parameterisations lie above the measurements. Other parameterisations of the pion structure 
function H44U45I were ruled out by previous HI and ZEUS measurements JHEl as they show a 
much flatter behaviour as a function of (3. 

The comparison of the measured F 2 N /T n and the parameterisations of the pion structure 
function is affected by uncertainties on the pion flux normalisation, as explained above. It may 
also depend on absorptive corrections [|46]449ll , which are not taken into account in this analysis. 
Neutron absorption may occur through rescattering which transforms the neutron into a charged 
baryon or shifts the neutron to lower energy or higher p T . 

The results presented here are consistent with the previous measurement by the HI Collabo- 
ration [1J. A similar analysis has been published by the ZEUS Collaboration [2J. There is good 
agreement between the two cross section measurements. For the extraction of the pion structure 
function, the flux factor normalisations used in the ZEUS analysis are different from the one 
used here. Within the uncertainties of the normalisation the HI and ZEUS results agree. 



5 Summary 

The cross section for leading neutron production in deep-inelastic positron-proton scattering 
da/dxL and the semi-inclusive structure function F 2 N (Q 2 , x,xl) are measured in the kine- 
matic region 6 < Q 2 < 100 GeV 2 , 1.5 • 10" 4 < x < 3 • 10~ 2 , 0.32 < x L < 0.95 and 
Pt < 0.2 GeV. The present measurements have experimental uncertainties of 10 to 15%. 

The measurements are well described by a Monte Carlo simulation including neutron pro- 
duction in fragmentation and neutrons produced from n + exchange, as predicted by the DJANGO 
and RAPGAP programs respectively. At large xl ~ 0.7 the 7r + -exchange process dominates. 

Within the measured kinematic range, the semi-inclusive structure function F 2 N and the 
inclusive structure function F 2 have similar (Q 2 ,x) behaviour, which is consistent with the 
hypothesis of limiting fragmentation. The dependence of F 2 LA ^ 3) on the variable is similar 
for all xl bins, in accordance with the expectation from proton vertex factorisation. The scaling 
violations observed in F 2 LN are similar in size and shape to those seen in the parameterisations 
of the inclusive structure functions of the pion and the proton. The data are used to estimate the 
structure function of the pion, up to uncertainties on the background contribution and the overall 
normalisation, in the framework of the one pion exchange model for the neutron kinematic range 
0.68 <x L < 0.77 and p T < 0.2 GeV. 
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Table 1: Bins in Q 2 , x and x L as used for the measurement of the semi-inclusive structure 
function F 2 iJV ^. 



x L range 


da/dxL [nb] 




da/dxL [nb] 






p T <x L - 0.69 GeV 


p T < 0.2 GeV 


0.30 


^0.37 


10.6 ± 0.11 


± 


1.3 


8.21 


± 0.09 ± 


0.99 


0.37 


^0.44 


12.0 ± 0.11 


± 


1.1 


6.75 


± 0.07 ± 


0.53 


0.44 


-r 0.51 


14.0 ± 0.11 


± 


1.3 


6.26 


± 0.06 ± 


0.45 


0.51 


^0.58 


15.6 ± 0.11 


± 


1.4 


5.92 


± 0.05 ± 


0.46 


0.58 


^0.65 


16.5 ± 0.11 


± 


1.4 


5.81 


± 0.05 ± 


0.42 


0.65 


^0.72 


16.7 ± 0.10 


± 


1.4 


5.71 


± 0.05 ± 


0.40 


0.72 


^0.79 


16.0 ± 0.10 


± 


1.3 


5.57 


± 0.05 ± 


0.31 


0.79 


^0.86 


12.4 ± 0.08 


± 


1.2 


4.54 


± 0.04 ± 


0.25 


0.86 


^0.93 


7.3 ± 0.06 


± 


1.2 


2.57 


± 0.03 ± 


0.38 


0.93 


-=- 1.00 


2.3 ± 0.03 


± 


0.9 


0.63 


± 0.01 ± 


0.26 



Table 2: Differential cross section da/dxL of leading neutron production in deep-inelastic 
scattering in the kinematic range 6 < Q 2 < 100 GeV 2 , 1.5 • 10~ 4 < x < 3 • 10~ 2 and 
0.32 < x L < 0.95. The first uncertainty is statistical and the second systematic. Normalisation 
uncertainties of 5% are not included. 
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Table 3: The semi-inclusive structure function F 2 (Q 2 ,x,xl), for neutrons with 
Pt < 0.2 GeV. The bin centre values of Q 2 , x, x L and (5 in the corresponding bins are also 
given. The first uncertainty is statistical and the second systematic. Normalisation uncertainties 
of 5% are not included. 
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Q 2 [GeV 2 ] x L p F 2 LJV(3) 
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Table 4: The semi-inclusive structure function F 2 ^'(Q 2 ,x,xl), for neutrons with 
Pt < 0.2 GeV, continued from Tabled 
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Q 2 [GeV 2 ] x L p F 2 LJV(3) 
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Table 5: The semi-inclusive structure function F 2 ^'(Q 2 ,x,xl), for neutrons with 
Pt < 0.2 GeV, continued from Table 
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Q 2 [GeV 2 ] x L p F 2 LJV(3) 



37 


9 


87 


10" 


3 





365 


1 


55 


10" 


2 





0469 


± 





0057 


± 





0046 


37 


9 


87 


10" 


3 





455 


1 


81 


10" 


2 





0446 


± 





0046 


± 





0030 


37 


9 


87 


10" 


3 





545 


2 


17 


10" 


2 





0390 


± 





0037 


± 





0038 


37 


9 


87 


10" 


3 





635 


2 


70 


10" 


2 





0372 


± 





0032 


± 





0023 


37 


9 


87 


10" 


3 





725 


3 


59 


10" 


2 





0335 


± 





0028 


± 





0019 


37 


9 


87 


10" 


3 





815 


5 


34 


10" 


2 





0326 


± 





0028 


± 





0018 


37 


9 


87 


10" 


3 





905 


1 


04 


10" 


1 





0138 


± 





0017 


± 





0023 


55 


2 


17 


10" 


3 





365 


3 


42 


10" 


3 





0962 


± 





0125 


± 





0074 


55 


2 


17 


10" 


3 





455 


3 


99 


10" 


3 





0621 


± 





0072 


± 





0074 


55 


2 


17 


10" 


3 





545 


4 


77 


10" 


3 





0727 


± 





0072 


± 





0045 


55 


2 


17 


10" 


3 





635 


5 


95 


10" 


3 





0622 


± 





0057 


± 





0034 


55 


2 


17 


10" 


3 





725 


7 


90 


10" 


3 





0540 


± 





0048 


± 





0041 


55 


2 


17 


10" 


3 





815 


1 


17 


10" 


2 





0486 


± 





0048 


± 





0027 


55 


2 


17 


10" 


3 





905 


2 


29 


10" 


2 





0204 


± 





0029 


± 





0039 


55 


4 


63 


10" 


3 





365 


7 


29 


10" 


3 





0695 


± 





0093 


± 





0052 


55 


4 


63 


10" 


3 





455 


8 


50 


10" 


3 





0614 


± 





0069 


± 





0052 


55 


4 


63 


10" 


3 





545 


1 


02 


10" 


2 





0448 


± 





0048 


± 





0032 


55 


4 


63 


10" 


3 





635 


1 


27 


10" 


2 





0482 


± 





0048 


± 





0035 


55 


4 


63 


10" 


3 





725 


1 


68 


10" 


2 





0498 


± 





0046 


± 





0035 


55 


4 


63 


10" 


3 





815 


2 


50 


10" 


2 





0415 


± 





0042 


± 





0034 


55 


4 


63 


10" 


3 





905 


4 


87 


10" 


2 





0164 


± 





0023 


± 





0022 


55 


9 


87 


10" 


3 





365 


1 


55 


10" 


2 





0592 


± 





0084 


± 





0048 


55 


9 


87 


10" 


3 





455 


1 


81 


10" 


2 





0429 


± 





0053 


± 





0027 


55 


9 


87 


10" 


3 





545 


2 


17 


10" 


2 





0409 


± 





0046 


± 





0030 


55 


9 


87 


10" 


3 





635 


2 


70 


10" 


2 





0372 


± 





0038 


± 





0027 


55 


9 


87 


10" 


3 





725 


3 


59 


10" 


2 





0338 


± 





0034 


± 





0022 


55 


9 


87 


10" 


3 





815 


5 


34 


10" 


2 





0304 


± 





0032 


± 





0018 


55 


9 


87 


10" 


3 





905 


1 


04 


10" 


1 





0148 


± 





0022 


± 





0037 


55 


2 


10 


10" 


2 





365 


3 


31 


10" 


2 





0515 


± 





0089 


± 





0046 


55 


2 


10 


10" 


2 





455 


3 


86 


10" 


2 





0347 


± 





0051 


± 





0021 


55 


2 


10 


10" 


2 





545 


4 


62 


10" 


2 





0286 


± 





0038 


± 





0038 


55 


2 


10 


10" 


2 





635 


5 


76 


10" 


2 





0334 


± 





0040 


± 





0018 


55 


2 


10 


10" 


2 





725 


7 


65 


10" 


2 





0303 


± 





0036 


± 





0020 


55 


2 


10 


10" 


2 





815 


1 


14 


10" 


1 





0266 


± 





0033 


± 





0023 


55 


2 


10 


10" 


2 





905 


2 


21 


10" 


1 





0098 


± 





0018 


± 





0018 


82 


4 


63 


10" 


3 





365 


7 


29 


10" 


3 





0588 


± 





0099 


± 





0045 


82 


4 


63 


10" 


3 





455 


8 


50 


10" 


3 





0537 


± 





0080 


± 





0036 


82 


4 


63 


10" 


3 





545 


1 


02 


10" 


2 





0465 


± 





0061 


± 





0036 


82 


4 


63 


10" 


3 





635 


1 


27 


10" 


2 





0520 


± 





0063 


± 





0052 


82 


4 


63 


10" 


3 





725 


1 


68 


10" 


2 





0480 


± 





0056 


± 





0033 


82 


4 


63 


10" 


3 





815 


2 


50 


10" 


2 





0417 


± 





0053 


± 





0040 


82 


4 


63 


10" 


3 





905 


4 


87 


10" 


2 





0180 


± 





0032 


± 





0028 


82 


9 


87 


10" 


3 





365 


1 


55 


10" 


2 





0484 


± 





0085 


± 





0041 


82 


9 


87 


10" 


3 





455 


1 


81 


10" 


2 





0514 


± 





0077 


± 





0036 


82 


9 


87 


10" 


3 





545 


2 


17 


10" 


2 





0366 


± 





0052 


± 





0025 


82 


9 


87 


10" 


3 





635 


2 


70 


10" 


2 





0507 


± 





0064 


± 





0038 


82 


9 


87 


10" 


3 





725 


3 


59 


10" 


2 





0361 


± 





0044 


± 





0027 




y 


07 
O 1 




3 


n 
u 


8 1 K 


K 



04 


in- 
Ill 


2 


u 


Uol 1 


zt 


u 




_1_ 

zt 


n 
u 




82 


9 


87 


10" 


3 





905 


1 


04 


10" 


1 





0170 


± 





0030 


± 





0027 


82 


2 


10 


10" 


2 





365 


3 


31 


10" 


2 





0477 


± 





0093 


± 





0040 


82 


2 


10 


10" 


2 





455 


3 


86 


10" 


2 





0319 


± 





0052 


± 





0025 


82 


2 


10 


10" 


2 





545 


4 


62 


10" 


2 





0366 


± 





0053 


± 





0024 


82 


2 


10 


10" 


2 





635 


5 


76 


10" 


2 





0372 


± 





0050 


± 





0030 


82 


2 


10 


10" 


2 





725 


7 


65 


10" 


2 





0443 


± 





0057 


± 





0028 


82 


2 


10 


10" 


2 





815 


1 


14 


10" 


1 





0241 


± 





0035 


± 





0020 


82 


2 


10 


10" 


2 





905 


2 


21 


10" 


1 





0132 


± 





0026 


± 





0027 



Table 6: The semi-inclusive structure function F 2 (Q 2 ,x,xl), for neutrons with 
Pt < 0.2 GeV, continued from Tabled 
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Figure 3: The observed neutron energy (a) and transverse momentum (b) distributions in the 
kinematic range 6 < Q 2 < 100 GeV 2 and 1.5-10~ 4 < x < 3-10 -2 . The data are compared to the 
predictions of RAPGAP-7T (dashed line) and DJANGO (dotted line) Monte Carlo simulations. 
Also shown is a weighted combination of those two simulations (full line), as described in 
Sect.H3 
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Figure 4: The cross section as a function of the fractional energy of the neutron xl in the 
kinematic range 6 < Q 2 < 100 GeV 2 and 1.5 ■ 10~ 4 < x < 3 • 10~ 2 . The transverse momentum 
of the neutron is restricted to p T < xl ■ 0.69 GeV (a) and pr < 0.2 GeV (b). The data are 
compared to the predictions of RAPGAP-7T (dashed line), DJANGO (dotted line) and LEPTO- 
SCI-GAL (dash-dotted line) Monte Carlo simulations. Also shown is a weighted combination 
of RAPGAP-tt and DJANGO simulations (full line), as described in Sect. [3^1 
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Figure 5: The semi-inclusive structure function F% N (Q 2 , x, xl), for neutrons with 
p T < 0.2 GeV, compared to the predictions of RAPGAP-7T (dashed line) and DJANGO (dotted 
line) Monte Carlo simulations. Also shown is a weighted combination of those two simulations 
(full line), as described in Sect. I3.4I 
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Figure 6: The ratio of the semi-inclusive structure function (Q 2 , x, for neutrons with 

p T < 0.2 GeV, to the proton structure function F 2 (Q 2 , x) obtained from the H1PDF2009 fit to 
inclusive DIS data ll26ll . The lines show the average value for a given xl bin. 
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Figure 7: The semi-inclusive structure function F 5 f JV *- 3 ''(Q 2 , /3, x£), for neutrons with 
Pt < 0.2 GeV, shown as a function of (3 in bins of Q 2 and xl- The lines are the results of 
the fit with a function c(xl) ■ /3 _A( - Q as described in Sect. HI 
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Figure 8: The semi-inclusive structure function F 2 LJV ^ 3 ' ) , for neutrons with p T < 0.2 GeV, di- 
vided by the pion flux T n integrated over t at the central value x L = 0.73, shown as a function of 
Q 2 in bins of (3. The pion flux is defined in Eq.|7J The data are compared to two different param- 
eterisations of the pion structure function F% [|42ll43l and to the H1PDF2009 parameterisation 
of the proton structure function 11261 . which has been scaled by 2/3. 
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Figure 9: The semi-inclusive structure function F^^, for neutrons with px < 0.2 GeV, di- 
vided by the pion flux rV integrated over t at the central value xl = 0.73, shown as a function 
of (3 in bins of Q 2 . The pion flux is defined in Eq. [7] The data are compared to two different 
parameterisations of the pion structure function F% ||42ll43l and to the H1PDF2009 parameter- 
isation of the proton structure function 11261 . which has been scaled by 2/3. The contribution of 
neutrons from fragmentation, as predicted by DJANGO and scaled by a factor 1 .2, as described 
in Sect. 13 .41 is indicated. 
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